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CHANGE POINT TESTING FOR THE DRIFT PARAMETERS OF A 
PERIODIC MEAN REVERSION PROCESS 

HEROLD DEHLING, BRICE FRANKE, THOMAS KOTT, AND REG KULPERGER 

CS| . Abstract. In this paper wc investigate the problem of detecting a ehange in the drift 

parameters of a generahzed Ornstein-Uhlenbeek proeess whieh is defined as the solution of 

^ ■ dXt = (Lit) - aXt)dt + adBt 

^ ■ and which is observed in continuous time. We derive an explicit representation of the 

_ O , generalized likelihood ratio test statistic assuming that the mean reversion function L{t) 

■ is a finite linear combination of known basis functions. In the case of a periodic mean 

reversion function, we determine the asymptotic distribution of the test statistic under the 
null hypothesis. 



H 

c3. 1. Introduction 

The problem of testing for a change in the parameters of a stochastic process has been an 
important issue in statistical inference for a long time. Initially investigated for i.i.d. data, 
^ ■ change point analysis has more recently been extended to time series of dependent data. For 

O . a general review of change-point analysis, see e.g. the book by Csorgo and Horvath [3]. 

In the present paper, we investigate the problem of detecting changes in the parameters of 
O I a diffusion process. Diffusion processes are a popular and widely studied class of models with 

applications in economics, finance, physics and engineering. Statistical inference for diffusion 
processes has been investigated by many authors, see e.g. the monographs by Liptser and 
Shiryaev [12] and by Kutoyants [11]. However, change-point analysis for diffusion processes 
has found little attention up to now. 

In our paper, we focus on change-point analysis for a special class of diffusion processes, 
^ ! namely for so-called generalized Ornstein-Uhlenbeek processes. These processes are defined 

d I as solutions to the stochastic differential equation 

(1) dXt = {Lit) - aXt)dt + adBt, t > 0, 

where a and a are positive constants and where the mean- reversion function L{t) is non- 
random. {Bt)t>o denotes standard Brownian motion and Xq is a square-integrable real- valued 
random variable that is independent of {Bt)t>Q. If L{t) = /x is a constant, we obtain the 
classical Ornstein-Uhlenbeek process, introduced by Ornstein and Uhlenbeck [14] . Ornstein- 
Uhlenbeek processes are popular models for prices of commodities that exhibit a trend of 
reversion to a fixed mean level. Generalized Ornstein-Uhlenbeek processes can be used as 
models for the evolution of prices with a trend or seasonal component L{t). 

Dehling, Franke and Kott [5] have studied the problem of parameter estimation of a gener- 
alized Ornstein-Uhlenbeek process if the mean-reversion function L{t) is a linear combination 
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of known basis functions ipi{t), . . . , (pp{t), i.e. when 

p 

(2) L{t) = J2^^^f^it) 

i=l 

In this model, the unknown parameter vector is ^ = (/ii, . . . ,fip,ay. We denote the corre- 
sponding parameter space by 0, and observe that 

= RP X (0,oo). 

As is usual in the statistical inference for the drift of a time-continuously observed diffusion 
process, the diffusion parameter a is supposed to be known. This can be justified by the 
fact that the volatility a can be computed by the quadratic variation of the process. 

We are interested in testing whether there is a change in the values of the parameters 
fii, . . . , fip and a, in the time interval [0, T] during which the process is observed. In the first 
step, we will consider this problem assuming that the change-point r G (0, T) is known. For 
the asymptotic analysis, when T — oo, we write t = sT, where s G (0, 1) is known. The 
generalized Ornstein-Uhlenbeck process with change-point r = s T is given by 

(3) dXt = {S{e,t,Xt)l{t<r} + S{e\t,Xt)l{,^,y)dt + adB^, < t < T, 
where 

p 

(4) S{9,t,Xt) = J2f'^f^it)-aXt, 

i=l 

and where 1^ denotes the indicator function of the set A. The test problem of interest can 
be formulated as 

(5) Hq : 6 = 6' (no change) vs. Ha : 6 ^ 6' (change at time point r) . 

We want to study the generalized likelihood ratio test for this test problem. We denote 
by Px the measure induced by the observable realizations X'^ = {Xt, < t < T} on 
the measurable space {C[0,T],B[0,T]), C[0,T] being the space of continuous, real-valued 
functions on [0,T] and S[0,T] the associated Borel cx-field. Moreover, let Pb be the measure 
generated by the Brownian motion on (C[0, T], ;B[0, T]). Then the likelihood function C 
of observations X'^ of the process with stochastic differential (3) is defined as the Radon- 
Nikodym derivative, i.e. 

C{6,6',X^) ■.= ^iX^)- 
The generalized likelihood ratio TZ{X^) is given by 

sup.ge i^{O.O.X^) 



Note that this likelihood ratio depends on the suspected change point r = sT, where s G 
(0, 1). Eventually, we will study the log-transformed likelihood ratio 

Aris) :=-21og(7^(X^)). 

We will give an explicit expression of the process (At(s))o<s<i and study the asymptotic 
distribution of this process as T — oo. 

The outline of the paper is as follows. In Section 2, we will first derive an explicit repre- 
sentation of the log-transformed generalized likelihood ratio test statistic (At(s))o<s<i. We 
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then formulate two theorems concerning the asymptotic distributions of sup-norm function- 
al of (A7-(s))o<s<i- It turns out that suPq<j,<]^ Aj'(s) does not have a non-degenerate hmit 
distribution, as T — cxd. In Theorem 1, we will prove convergence of supg^<g<g2 Ar(s), when 
< Si < S2 < 1 are fixed constants. In Theorem 2, we will show that there exist centering 
and norming sequences qt and bx such that (supo<s<i At(s) — bT)/aT converges towards 
an extreme value distribution. The proofs of these theorems are given in Section 3 and 
Section 4, respectively. 



2. Generalized likelihood ratio test 

In this section, we will derive an explicit representation of At{s). In order to do so, we 
need to calculate the maxima in the numerator and denominator in (6). Note that this is 
achieved by the corresponding maximum likelihood estimators. A corollary to Girsanov's 
theorem, see Theorem 7.6 in Lipster and Shiryayev [12], gives an explicit expression of the 
likelihood function of a diffusion process provided that 

(7) ^^lo -^(^'^'^t)'^^ < °^ 
for all < T < oo and all 6. 

Lemma 2.1. Let C{9,X'^) denote the likelihood function of the observations X'^ = {Xt,0 < 
t < T} of the generalized Ornstein- Uhlenbeck process {Xt)t>o, defined in (1), with mean 
reversion function (2). If the drift term (4) satisfies condition (7) then 

aigmax C{6, X'^) = 6ml = Q^^Rt- 

9 

Here Qt G M(p+i)x(p+i) and Rt G M^+i are defined as 

- y -dr br 
Rt = <^i(t)dXt,...,^ ^,{t)dXt,- XtdX^j , 

where Gt = {Jq '^j{t)<^k{t)dt)i<j^k<p e ar = {J^ ipi{t)Xtdt, . . . , ipp{t)Xtdty and 

bT = J^X^dt. 

Remark 1. Note that the integrals in Rt can be rewritten as 

iPi{t)dXt = / V^{t){L{t) - aXt) dt + a ^i{t)dBt 
Jo Jo 

where the latter is a well-defined Ito integral. 

Proof. The likelihood function £ of a general diffusion process 

dXt = S{9, t, Xt)dt + adBt, < t < T, 

is given by 

(8) C{e,X^) = ^{X^)=exp(^ s{9,t,X,)dX,-^ S{e,t,Xtfdt\ 
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if condition (7) is fulfilled; see Theorem 7.6 in Lipster and Shiryayev [12]. The maximum 
likelihood estimator is defined as the maximum of the functional 9 ^ C{6, X^) and the 
partial derivatives of the logarithm of this functional are 

o -i p'Y' o -1 pT o 

(9) _in(£(^,x^)) = - / —sie,t,x^)dXt-- sie,t,Xt)^sie,t,x,)dt. 

The derivatives of the drift function specified in (4) can be computed to be 

Setting the partial derivatives of the log-likelihood function in (9) equal zero gives a system 
of linear equations which yields the assertion. □ 

Due to the linearity of the drift term, the log-likelihood function of the process (3) is given 

by 

in(£(^,^',x^)) = sie,t,Xt)dXt + sie',t,Xt)dx, 

Hence, defining X^'^ = {Xt, t <t<T}, we can write the generalized likelihood ratio (6) as 
^ ' ^ ^ snpg,C{9*,X-)snps,C{9',X-'T) 

where CiO^X'^) is given in (8) with drift function specified in (4). The terms C{6*,X'^) and 
C{6',X'^'^) are defined analogously as integrals with integration regions to r and r to T, 
respectively. It follows from Lemma 2.1 that 

(11) n{x^) = . ^^^''"^'^^^ 

c{ei,,,x-)c{e]:^,,x-n 

where the maximum likelihood estimates Oj^j^^, 6\jj^ and 9]j[i^ are computed from the total, 
the pre- and post-change sample, respectively. This representation of the likelihood ratio is 
used to prove the following proposition. 

Proposition 2.2. The log-transformed generalized likelihood ratio test statistic At(s) = 
— 2 ln(7^(X'^)) of the test problem (5) can he represented under the null hypothesis as 

Kt{s) = -R'tQt^Rt + KQr^Rr + {Rt ' RrYiQx ' Qr^^Rx - Rr) 

where Qt is given in Lemma 2.1 and 

i?T=(^^ Vi{t)dBt,..., iPp{t)dBt,- XtdBt^ . 

Proof. Our aim is to compute an explicit expression of the ratio given in (11). Note that the 
likelihood function in the numerator of (10) can be represented as 

L{9,X^) = exp (^O'Rt - ^^'QtO 
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where Rt and Qt are given in Lemma 2.1. Denoting by 6q the true value of 6, the represen- 
tations 

^ML ~ Qj^Rt and Rx = Qt^o ~I" ^Rti 
where the latter can be obtained by plugging in the initial SDE (1) and (2), lead to 

L(C^,X^) = exp (^^i?^Q^ii?r^ 

= exp {^^^^qQt^^o + —R^do + -^RtQt^ ■ 

The same procedure yields an analog expression for L{9\ji^,X^). The additivity of the 
integrals provides 

^ML = ^Qt - QrY^iR-T - Rt) 



and 

such that 



Rt — Rt = {Qt — Qt)0o + o"(-Rt — Rt 



L{e]^l,x^'^) = exp{—el{QT-QT)eo + -{RT-RT)% 



^^0l{QT-QT)eo + -. 

Za a 

+i(i?T - Rt)\Qt - Qt)-\Rt - Rt)^ ■ 

Under the null hypothesis, cancelation of several terms in (11) proves the assertion. □ 

For the rest of our investigations, we study periodic functions 
(12) ^,it + u) = ^,it) 

where u is the period observed in the data. Under the null hypothesis of no change, this 
implies periodicity of the mean reversion function, i.e. L{t + u) = L{t). We assume that 
we observe the process over some integer multiple of periods, i.e. T = nu, G N. By 
Gram-Schmidt orthogonalization we may assume without loss of generality that the basis 
functions ipi{t), . . . , (pp(t) form an orthonormal system in L'^{[0, u], -d\), i.e. that 



(13) / Vj{t)Mt)dt 



j = k 
0, J k. 

Under these assumptions, the matrix Qt appearing in the test statistic Ar(s), see Proposition 
2.2, simplifies to 

T Jpxp O-T 

atp Bt 



Q 



T 



Theorem 1. Let = {Xt,0 < t < T} be observations of the mean reversion process 
(1) with mean reversion function of the form (2), satisfying (12) and (13). Denote by 
Ax'(s) = —2\n{TZ{X^)) the log-transformed generalized likelihood ratio test statistic for the 
test problem (5). Then, for any fixed < Si < S2 < 1? under the null hypothesis, 

_^ \\w{s)-sw{i)r 

sup IyT\S) > sup r — 

se[si,s2] se[si,s2] '^v-'- ~ ^) 

as T ^ oo. Here \\ ■ \\ denotes the Euclidean norm and W is a {p+ 1)- dimensional standard 
Brownian motion. 
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The result stated in Theorem 1 is not satisfactory in apphcation: First, it is not clear 
how to choose the interval [si,S2] potentially containing a change point if no information 
about the location of the change is available. Second, the distribution of the limit which is 
the squared length of a mult i- dimensional Brownian bridge is not explicitly given such that 
further analysis or a simulation study are necessary in order to specify the limit distribution 
explicitly. In order to avoid this inconvenience it is possible to consider the exact test statistic 
supo<s<i Ar(s). It turns out that, by means of some appropriate normalizing terms ax and 
Bt, the distribution of the expression (supQ<g<]^ At(s) — hT)/aT converges to the Gumbel 
distribution. 

Theorem 2. Under the same assumptions as in Theorem 1 it holds under the null hypothesis 
that 

( sup At(s) — hr) /ar — > G, 

0<s<l 

as T oo, where G denotes a real-valued random variable satisfying 

P(G <x)= exp(-2e-^/2)_ 



Herebr = (2 Inln ^ + £±i Inlnln ^ - In r(2±i)) V(2 Inln ar = ^Jbr/ {2\n\n^) where T 
is the gamma function. 

3. Proof of Theorem 1 

Before we can complete the proof of Theorem 1, we have to establish some auxiliary results. 
First, we will study the asymptotic behavior of At{s) in the case of a periodic mean reversion 
function, see (12) and (13). Note that by Proposition 2.2 we have the representation 

At(s) = —R^rpQrp^ Rt + I^grpQsT^sT ^ {Rt " RsxYiQr — Qst) ^{Rt — Rst)- 
The first term, 

RtQt^Rt = -^Rt{i^Qt^ v^^^' 

has already been studied by Dehling, Franke and Kott [5]. The following proposition sum- 
marizes the results of Proposition 4.5, 5.1 and 5.2 in Dehling et al. [5]. 



Proposition 3.1. We have that 



and 



Rt N{0, S) 



—Qt almost surely, 



as T ^ oo. The matrix S is given by 
(14) S 



A* CO 



where Aj = Lpi{t)h{t)dt, i = 1, . . . ,p, co = J^^ {h{t)ydt + ^ and where h : [0, oo) — M zs 
defined by 

(15) /i(t) = e-"*^/ij- / e">j(s)ds. 

j=l J-oo 
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Here, N{0, S) denotes a normally distributed random vector with zero-mean and covariance 
matrix S. 

Now we want to investigate the second term of At{s) which we rewrite as 

1 



^ItQsT^sT 



We will show that the process (^^Rst) ^^^^^ converges in distribution to a Gaussian process 

on [si, S2], and that ^Qst converges in probability uniformly on [si, S2]. 

We need the following functional version of the asymptotic normality proved in [5] 

Proposition 3.2. As T ^ 00, the sequence of processes {R^^)o<t<i, where 

1 



T 



R 



tT, 



converges in distribution to a {p + 1)- dimensional Wiener-process R* with R* ~ A^(0, sS) 
and where S is defined in (14)- Thus the covariance function of R* is of the form 

Cov {R*{s),R*{t)) = (s At)S,,; fort,j 
Proof. Note that the vector-valued processes 



R 



in) 



nut 



ipi{s)dBs 



nut 



iPp{s)dBs 



nut 



Xi^dBf, 



are martingales with respect to the filtrations 



(n) 



<^{Bs',s < nut). The associated 



covariance processes are given by 

)t — I 1 rnut J 1 rnut , i rnut -ia j 

As was shown in [5] (see p. 184), these matrices converge for n — 00 almost surely towards 
the matrix 

^ t (p lipids ... t (fitppds t Jq ipihds ^ 



V 



2a 



J 



t ^ ippipids ... t Jl^ ippippds t ipphds 
t Jl^ (pihds ... t (pphds t (^f^ Ji^ds -\- 

The functional central limit theorem for continuous martingales (p. 339 in [6]) now implies 
that the sequence of continuous J-'j"^-martingales i?^"^ converges in distribution toward the 
unique continuous Gaussian martingale with covariance function tE. □ 

Proposition 3.3. Let Qt be defined as in Lemma 2.1. Then, as T ^ 00, 

1 



T 



Q 



sT 



almost surely uniformly on [0, 1] , where S is given in (14). 
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Proof. By Proposition 3.1, we know that, almost surely, ^Qt — )■ S as T — )■ oo . Thus, given 
e > 0, there exists a Tq such that for all T >Tq 

II^Qt-S|| <e. 

Let B := supo<s<T„ \\Qs\\- Then we get for any T > Tq and Tq/T < s < 1 



For s < To/T we obtain 



\j;Qts - sS|| = s\\^Qts - S|| < e. 



\^QTs-sm<^B + ^\\j:\\<e, 



for T large enough. The last two inequalities together show that for T large enough, we have 
II ^^Qts ~ -5^11 < e, and this proves the statement of the proposition. □ 

We can finally finish the proof of Theorem 1. 

Proof of Theorem 1. By Slutsky's theorem and Propositions 3.1, 3.2 and 3.3 we obtain 

At('S) = —B!^Q7p^RT + R\rpQ~^RsT + {RT — RsTy{QT — QsT)~^{.RT — RsT) 

^ -\\W{1)\\' 



\\W{s)r , \\W{l)-W{s)r 



s 1 — s 

\\w{s)-sw{i)r 

s(l-s) 

in C[si,S2]. Here we have used the fact that the process W{t) := is a Brownian 

motion with covariance matrix Jp+i, where /p+i is the (p+ 1) x 1) identity matrix, and 
that (-RJ')*S~^i?j = ||M^(t)||^. The assertion about the supremum of Kt{s) is justified by the 
continuous mapping theorem. □ 

4. Proof of Theorem 2 

The proof of Theorem 2 is motivated by the proof of an analogous result for discrete time 
AR processes, given by Davis et al. [4]. We need the following result which is proved in in 
Davis et al. [4] and which relies on Lemma 2.2 in Horvath [7]). 

Proposition 4.1 (Corollary A. 2 in Davis et al. [4]). Let Yi,Y2, . . . be an i.i.d. sequence 
of {p + 1) - dimensional random vectors with lE[Yi] = and ]E[YiY-l^] = Ip+i- Define S, 

Si=i ^i- If ^^^i<i<p+i lE|^i,iP''''' < oo for some r > 0, then 

( max \\Sk\\'^ - bn)/an G* , 

^ l<k<n ' ' ^ 

as n ^ oo, where G* denotes a real-valued random variable satisfying 

P{G* <x)= exp(-e-"/'). 



k 



Thereby, it is bn = (2 Inlnn + 2±i Inln Inn - lnr(2±i)) /(21nlnn), a„ = ^/bn/{2\n\nn). 
Recall that 

At(s) = —R^rpQ^^Rrp + R^rpQ^^Rgrp -|- (^Rrp — RsT^'^Qt " Q ^ {Rt " RsT) 

and T = raz/, v fixed. Let us assume for a moment that u = 1. We write A„(s), -R„ and Qr 
for At(s), Rt and Qt, respectively, since the asymptotic framework is n — )■ oo. 
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Assume further that the solution of the stochastic differential equation (1) is stationnary. 
Note that this is not a constraint for our purposes since we may alternatively consider the 
stationary process 

Xt := hit) + Zt 



with 



and 



where 



ft 

-at 



h{t) := e""* / e'''h{s)ds 



ft 



oo 



is bilateral Brownian motion defined through two independent standard Brownian motions 
{Bt)t>o and {Bt)t>o- The process {Xt)t>o is stationary and does not depend on Xq. Fur- 
thermore, it was proved in [5] (Lemma 4.4) that one has \Xt — Xt\ — 0, almost surely, as 
t — J- oo. 

The following proposition is essential for the proof of Theorem 2. First, define for two 
cr-algebras A and B the quantities 



and 



a{A,B)= sup \P{AnB) -P{A)P{B)\ 

A£A,B£B 



p{A,B)= sup Corr(F,G). 

F6L2(^,p),GeL2(B,P) 



It is known that 

a{A,B)<^p{A,B). 

For a stationnary sequence of random variables {Ck)ken define the mixing coefficient by 

af(n) = sup a(a{Q; i < k), cr{Q; i > k + n)). 

fceN 

The sequence (Cfc)fceN is called strongly mixing if — as n — )■ oo. 

Proposition 4.2. The sequence of random vectors {rk)keN defined by 

(pk pk j'k \ * 

/ ipiit)dBt, . . . , ipp{t)dBt,- XtdBt) 
Jk~l Jk-l Jk~l J 

is strongly mixing with mixing coefficient a of order 

Proof. Define the C[0, l]-valued stochastic process (^^^■')fceN by 

t{k) ._ ( -^^^'^ \ ■- ( {^t+k)t(i[Q,l\ 

^ V B^'^ )'~\ iBt+k)tem 
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The process {^^^^)k£N is both a Markov and a Gaussian process. Hence, by making use of 
the Markov property and by applying the correlation inequality for Gaussian processes from 
Lemma 4.3 we obtain 

a^{n) < supa(o-(^^''); A; < m),cr(^(''); A; > m + ra)) 

teM 

= supa(cT(e("^)),(T(e^"^+"^)) 

< sup sup Corr((a,5('")) + (6,X(™)), (c,5(™+")) + 

meNa,6,c,de-L2[o,l] 

< e-"(""^) sup Corr((a,5W) + (6,xW),(c,5(2)) + 

a,b,c,d,eL'^[0,l] 

where the last equality is stated in Lemma 4.4. Note that each ARk may be represented as 

where fk : C[0, 1] x C[0, 1] — R is a measurable function. Since the cx-algebra generated by 
fkiC^''^) : — )■ M is smaller or equal the cx-algebra generated by C,^''^ : Q — )■ C[0, 1] x C[0, 1], 
the bound for a^{n) established above is also valid for a/\ii{n). □ 

Lemma 4.3. Let [H, (., .)) be separable Hilbert-space and {X,Y) be a pair of random vari- 
ables with Gaussian joint law. Then one has 

p{a{X),a{Y)) < max Corr((a, X), (6, F)). 

a,b£H 

Proof. Let (ej)jgN be a system of orthonormal basis vectors for the Hilbert space H. If we 
set Vi := (X, Cj) and Wj := {Y,ej) then we have the representations 

oo oo 

X = J2 ViCi and Y = J2 ^fir 
i=i j=i 

Note that a{X) = a{Vi,V2,...) and a{Y) = a{Wi,W2, ...)■ It follows from Prop. 3.18 and 
Thm. 9.2 in [2] that 

p{a{X),a{Y)) = \un p{a{V^,..,Vn),aiWi,...,Wn)) 

n— >oo 

n n 

= lim sup Corr( N aiVi,} bjWj). 

"^°°ai,...,a„,bi,...,b„6M ^ 

Since the correlation is homogeneous we can assume without loss of generality that Yl^l ~ ^ 
and ^ = 1 holds. From this then follows 

oo oo 

p(cr(X),a(y)) < sup Corr(^aiV^,,^6,W^, 

(ai)!GN,(6j)jeN:Eaf=E''?=l i=l j=l 

< sup Con f {a, X),{b,Y) 

a,feeH:||a|| = ||6||=l ^ 

The second inequality follows since one has ioi a E H and := (a, Cj) that 

oo 

{a,X) = Y,(^^V^■ 
1=1 
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This finishes the proof of the lemma. □ 
Lemma 4.4. For all a, b,c,dE L'^[0, 1] we have 

Corr((a,5(™)) + (c,fi('"+")) + 

= e-"("-^)Corr((a,5('")) + (c,5(™+i)) + 

Proof. Since is independent from cr(-B''™'', X*^™^) we have that 

Corr((a,5('")) + (&,X(™)), (c,5(™+'^)) + (rf,X('"+"))) 
= Corr((a, E^™)) + {b, X^™)), (rf, X(™+"))) 

Cov((a,fi('")) + (^X^'")), (d,X('"+"))) 
v/Var((a,fi(")) + (6, XH)) v/Var((rf, 

Cov((a, + {b, X(™)), (rf, X(™+"))) 

v/Var((a,5(")) + (6, X^)) v/Var((c/, XW)) ' 

Note that we used the fact that the sequence (X^™'*)meN is stationary. 

We will use the fact that {Xt)t>rn is the unique solution of the SDE (1) with initial condition 
Xm to see that X*^™"*"") has the representation 



/•n+s 

X^+n+s = e-"("+^)X„ + h{n + s) + ae-"("+^) / e"^c^S^+,. 

Jo 



This representation follows from Lemma 4.2 in Dehling, Franke, Kott (2010). We use this 
fact to compute the covariance in the above formula: 



Cov((a,B(™)) + (6,X(™)), (d,X('^+"))) 



Gov (a,5(™)) + (6,X(™)), / c{s)e~'''Xmds 







+Cov| (a,5(^")) + (6,X(™)), / c(s)/i(n + s)c/s 







+e-""Cov|^(a,B(™)) + (6,X('")),^ c{s)ae-''' e""' dB^n+rds 
+e-"'^Cov ^(a, B^'")) + (6, X^'")), c(s)ae-"^ ^ 



6 dBr^iA^j-ds 



Note that the second term on the right vanishes, since the right entry in the covariance 
is deterministic. Further, the fourth term vanishes, since the Brownian increments on the 
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interval [m + 1, m + n + s] are independent with respect to a{B^'^\ X^"^^). We thus have 



e-""Cov ^(a, + {b, ^ 



c{s)e "'^Xmds 



+e-""Cov I (a, E^*")) + (6, X^""^), j c{s)ae' 



The result follows since we can do the same reasoning for n = 1. 



□ 



Corollary 4.5. There exists an iid-sequence of W'^^ -valued Gaussian random variables 
Ci; ^ G N such that for Uk ■= Yl'i=i d one has 

- Uk = 0{k^^'^'^) for some X > as k ^ oc. 

Proof. This follows from Proposition 4.2 and the theorem from Kuelbs and Philipp on strong 
approximation of mixing random sequences (see [9]). □ 

Remark 2. In the following we will denote by Fp+i the covariance matrix of the Gaussian 
random variable Ci- It then follows that the sequence of random variables Tp^-j^R[nt] / t E 
[0, 1] converges in distribution toward a. p+ 1-dimensional Brownian motion with covariance 
matric Ip+i- Here /p+i denotes the identity matrix with p+1 rows. It follows from Proposition 
3.2 that Fp+i = E. 

Proposition 4.6. For all 6 > one has as u ^ Q: 



lim sup P 

T-!>oo 



sup At(s) 

0<S<M 



sup E^grpQ ^^Rgrp 

Q<s<u 



and 



limsupPf sup At(s) - sup {Rt - RstY{Qt - Qst)~\Rt - Rst] 

T-s>oo V l-u<s<l l-u<s<l 



> arS -> 



Proof. It holds that 



(16) 



a^^l sup Aj'(s) — sup RItQsT^st] 

0<s<u 0<s<u 

< sup aj^^\AT{s) - RIj^QJ^^Rst] 

0<s<u 

= sup a^^\{RT — RsrYiQT — QsTy^iRx — Rst) — RtQt^Rt\ 

0<s<u 

ii'^w-«'W£_iwi)f 



sup 

0<s<u 



1 - s 

—J- 0, almost surely 



(as T — oo) 



CHANGE POINT TESTING FOR A PERIODIC MEAN REVERSION PROCESS 



13 



as n — 0, where the convergence in (16) is imphcated by the proof of Theorem 1 and the 
fact that — 1- Analogously, one has 



sup Ar(s)- sup {RT-RsT)\QT-QsTy\RT-RsT] 

l-U<S<l 1-U<S<1 

< sup cij} \ R!^gj'Q gj'RsT — Rlj^Qrp^RT^ 

l~u<s<l 



V 



sup 

1 — M<S<1 



as T — 7- ool 



— )■ 0, almost surely, 



as M — 7- 0. 



□ 



Proposition 4.7. Under the framework of Theorem 1 it holds under the null hypothesis that 

— ( sup R^tQ'tRsT - hr^ A G* 

OT \Q<s<u J 

and 

— ( sup {Rt-Rst)\Qt-Qst)-\Rt-Rst)] 

CtT \l-u<s<l J 

as n ^ oo, where G* denotes a real-valued random variable satisfying 

]P(G* <x)= exp(-e-"/2) 
and where ar and hx are given in Proposition 4-i- 

Proof. The reasoning follows the lines of the proof of remark A3 presented in [4] (see page 
297). We first note that 

RIt^pIiRsT - Uf^T]^~l^UlsT] = RIt^~Ii{RsT " U[sT]) + (^^r - Ul^T])'^pllU[sT]- 

The law of iterated logarithm implies t/j'^^jF"^^ = 0(([sT] log[sT])i/2) and Corollary 4.5 then 
implies -R*^r~^^ = 0(([sT] log[sT])^/^). Using those facts and Corollary 4.5 again yields 

(17) R'sT^^liRsT - UIt]^;1,U[st] = OiisT]'^'') 

for some A' > as T — )■ oo. 

Since by Proposition 3.3 one has Qst/sT — )■ Fp+i it follows that 



RItQstRsT TfRsT^pliRsT 



R 



^ ^TD'^fr Q^T\^^i RsT 

jj^Sl P+1 - ) i p+1 



This relation together with Equation (17) implies that as T — )■ oo one has 
(18) RItQJtRsT - Ul^f;l,Ui,T]/[sT] 0. 

Proposition 3.3 and the continuous mapping theorem yield 

V 



sT) 



1/2 



0. 



It thus follows that 



sup rCgT^sT^sT > sup . 

se[u,i] se[u,i] -5 



sup R\tQ~stRsT = Op(l) 
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Moreover we have 

sup Rgj'Q^rpRsT — )■ oo. 

s6(0,u] 

Thus with propabihty closer and closer to one the supremum is achieved in the interval (0, u] 
and not in [u, 1]. It then follows that 



(19) 



P( sup Rl^QjR^T= sup RItQsT^st) — > 1- 



We also note that for a fixed M > one has 
(20) 

and 
(21) 

Let 



P( sup RItQJRst= sup RItQJRst) — ^ 1- 
^se{o,i] se(A//T,i] ^ 

P( sup Ul^fpl^UlsT] = sup Ul^fpl^UlsT]) > 1- 

^sG(0,l] sG(A//T,l] ^ 



Rt{M):= sup Rl^Qj^Rl^- sup f/f,r]r;^if^[.T]/['^^]- 

se{M/T,i] se(Af/r,i] 



From Equation (18) we have 



(22) 



\Rt{M)\ 



sup RItQJ^RIt- sup Ulj^f;l,Ulj.J[sT] 

s6(M/n,l] se(Af/n,l] 



< sup 

sG(A//n,l] 



which goes to zero as M — > oo uniformly in T > u. It now follows from Eq.(19) and 
Proposition 4.1 that 



T-5. 



M 



lim W ( sup R\j.Qg^RsT < ctTa^ + &t 
lim P( sup R^rpQ'^Rgrp < + 6j. 

lim lim P( sup RItQJt-^sT ^ cltx + bT] 

i^oo T^oo V se(A//T,l] ^ 

lim hm P f sup U^T] KliUfsT] /[sT] < arx + br + Rt{M)) 



■ s£{M/T,l] 

hm P( sup ^f,T]r;^it^[.T]/[^2-] < arx + br 



= lim Pf sup \\S[sT]\\V[sT] < arx + br) — ^ exp ( - e"^'/^). 

This proves the first statement of the proposition. The second one is proved in an analogous 
way. □ 

Proof of Theorem 2. Since for fixed x G M one has a^x + 6^ — t- oo as T — oo it follows from 
Theorem 1 for all u G (0, 1/2) that 



p[ sup At{s) < arx + br) — > 1. 
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Therefore, one has as T — t- oo that 

P( sup Ar(s) < axx + fer) = f ( sup At{s) < a^x + 6^, sup A'r(s) < a^a; + Bt 

^0<s<l ' Vo<s<M \-u<s<\ 

By Proposition 4.6 this has for T — t- oo the same hmit as 

P sup RItQsT^sT < cltx + br, sup {Rt - Rst)\Qt - QstY^^Rt - Rst) < otx + 6t I • 

Proposition 4.7 yields that this last expression converges toward exp(— 2e~^'/^) since the two 
sequences 

sup Rl'pQgr^RsT 

0<s<u 

and 

sup (RT-RsTYiQT-QsTr'iRT-RsT) 

1-M<S<1 

are asymptotically independent by Proposition 4.2. □ 
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